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INTRODUCTION 


In  the  past  ten  years,  many  reports  have  been  published  concerning  the 
application  of  computer  programming  for  the  calibration  of  photographic  emulsions. 
Calculating  boards3  and  rules2  have  essentially  been  replaced  by  readouts  from 
large  and  small  computers.  Authors  have  generally  chosen  to  use  programs  that 
follow  graphical  methods  of  calibration, 3-8  although  some,  like  Margoshes  and 
Rasberry,8  take  a different  approach  by  assuming  a linear  relation  between  func- 
tions of  the  relative  exposure  and  of  the  microphotometer  readings.  Most  graphic 
methods  studied  appeared  to  be  unduly  complex;  for  example,  data  was  generally 
converted  to  Seidel  or  "modified"  Seidel  transformations,  the  use  of  the  prin- 
ciples involved  in  Kaiser  transformations  were  employed,  or  dedicated  computers 
were  used. 

Our  laboratory  has  striven  for  a program  that  would  follow  our  present, 
slightly  modified  ASTM  method  of  calibration10  which  is  simple  in  principle.  We 
sought  a computer  method  that  would  plot  preliminary  and  emulsion  calibration 
curves  so  that  a convenient  check  could  be  made  on  the  original  data  at  any  time 
desired.  By  retaining  the  ability  to  have  curves  printed  by  the  computer  one 
can  decide,  for  instance,  whether  more  points  are  needed  to  establish  the  prelim- 
inary curve,  as  it  must  be  remembered  that  the  accuracy  of  the  final  readout  of 
any  computer  program  is  limited  by  the  accuracy,  number,  and  selectiveness  of  the 
original  dark/ light  data  read  from  the  microphotometer. 

Using  our  program,  dark/light  data  are  inputted  to  obtain  the  emulsion 
calibration  equation.  This  equation  is  then  stored  either  in  the  memory  of  the 
calculator  or  on  recording  tape.  To  find  a relative  intensity,  percent  trans- 
mission readings  (one  from  the  element  in  question  and  the  other  from  an  internal 
standard)  are  entered  in  the  computer.  Troubleshooting  the  emulsion  calibration 
is  simplified  by  printing  the  preliminary  curve  with  the  original  data  super- 
imposed on  it  along  with  the  emulsion  calibration  curve. 

Two  similar  programs  were  originated  to  follow  our  existing  "hand  emulsion 
calibration"  method.  One  program  uses  time-sharing  computer  facilities  and  has 
the  capability  of  printing  charts  or  tables  of  percent  transmission  and  corre- 
sponding relative  intensities.  The  other  program  uses  a Hewlett-Packard  (H.P.) 
9830-A  calculator  with  an  8K  memory  storage.  An  X-Y  plotter  is  attached  to  plot 
curves  originated  by  the  data  from  the  computer. 

1.  KERBYSON,  J.  D.  A Calculating  Board  Incorporating  a Linear  Density  Transform . Spectrochimica  Acta,  v.  19,  1963,  p.  1335-1341. 

2.  CHURCHILL,  J.  R.  Techniques  of  Quantitative  Spectrograpliic  Analysis.  Ind.  Eng.  Chem.  Anal.  Ed.,  v.  16,  1944,  p.  653-670 
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in  Quantitative  Spectrometry.  Analytical  Chemistry,  v.  47.  March  1975,  p.  441-447. 
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Spectroscopy,  v.  10,  1965,  p.  33-37,  43. 
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in  Photographic  Microphotometry.  Spectrochimica  Acta,  v.  24H.  1969,  p.  497-513. 

10.  Methods  for  Emission  Spectrochcmical  Analysis.  American  Society  for  Testing  and  Materials.  1971. 


I 


PROCEDURE  FOR  EMULSION  PLATE  CALIBRATION 


Cassette,  rotation,  and  position  settings  on  the  spectrograph  were  adjusted 
so  that  the  wavelength  region  2750  to  3450  X could  be  photographed.  A solution 
containing  50  mg/ml  iron  was  sparked  for  a twenty-  and  forty-second  exposure. 
Details  concerning  the  parameters  used  for  excitation  and  conditions  set  on  the 
spectrograph  are  found  in  Appendix  A.  A seven-step  filter,  with  segments  in  the 
ratio  of  1:2,  was  placed  between  the  spark  and  the  entrance  slit  of  the  snectro- 
graph  to  stepwise  attenuate  the  exposure. 

Using  a densitometer,  lines  on  the  photographic  plate  in  the  wavelength 
region  3000  to  3200  A were  examined.  Transmittances  between  1.5%  and  92%  in  the 
66%  and  33%  portion  of  a line  (the  second  and  third  step  of  the  filter)  were 
listed  in  a table  labeled  DARK  percent  T and  LIGHT  percent  T.  The  preliminary 
curve  (PC)  is  generated  from  data  which  consists  of  about  50  line  pairs  plotted 
and  connected  on  linear  graph  paper.  Paper  15"  x 10"  in  size  was  labeled  DARK 
percent  T on  the  ordinate  and  LIGHT  percent  T on  the  abscissa.  Formation  of  this 
curve  is  the  first  stage  in  the  ASTM  method. 

The  second  stage  of  the  method  is  the  construction  of  the  emulsion  calibra- 
tion or  H and  D curve  and  it  is  the  equation  of  this  curve  that  is  ultimately 
used  in  relating  percent  transmission  to  relative  intensity.  To  construct  this 
curve  a new  table  was  established  using  the  two  parameters,  percent  T and  NUMBER. 
A DARK  percent  T point  was  chosen  on  the  preliminary  curve  as  low  in  value  as 
possible  but  above  the  intercept  of  the  curve  with  the  ordinate.  A corresponding 
LIGHT  percent  T point  was  obtained  and  this  value  was  written  under  percent  T 
and  7 under  NUMBER  as  in  the  table.  Allowing  the  value  of  the  LIGHT  percent  T 
to  equal  the  numerical  value  of  a DARK  percent  T point,  a new  LIGHT  percent  T 
value  was  obtained  from  the  PC.  In  the  table,  6 was  placed  under  NUMBER  and 
the  new  LIGHT  percent  T value  recorded.  This  method  was  continued  until  LIGHT 
percent  T values  were  obtained  for  the  corresponding  numbers  5,  4,  3,  and  2. 

Using  standard  size  8-1/2-  x 11-inch,  two  cvcle  semilog  paper,  the  ordinate 
or  log  scale  was  labeled  percent  T (1-100)  and  the  linear  portion  or  abscissa 
was  marked  off  from  left  to  right,  one  to  eight  with  ten  divisions  between  marks. 
This  axis  was  labeled  NUMBER.  Values  from  the  table  were  plotted  and  a straight 
line  was  drawn  connecting  the  two  points  corresponding  to  the  values  of  6 and  7 
and  extending  to  the  abscissa  axis.  After  establishment  of  this  line,  inter- 
mediary points  between  NUMBERS  are  taken  so  that  the  calibration  curve  can  be 
drawn  more  accurately.  This  was  accomplished  by  choosing  an  abscissa  value 
between  7 and  8,  for  example,  7.5,  and  obtaining  the  corresponding  percent  T 
value  from  the  straight  line  drawn  on  the  emulsion  calibration  curve.  This  value 
was  set  equal  to  a DARK  percent  T point  on  the  PC  and  referring  back  to  that 
curve  a corresponding  LIGHT  percent  T value  was  determined.  In  the  previously 
established  table,  place  the  value  6.5  under  NUMBER  (6.5  chosen  because  it  is 
one  less  than  7.5,  the  starting  value)  and  the  newly  obtained  LIGHT  percent  T 
value  under  percent  T.  Continuing,  let  this  new  LIGHT  percent  T value  equal  that 
of  a DARK  percent  T point  on  the  PC  and  obtain  a new  LIGHT  percent  T number 
from  ihe  PC.  In  the  table,  place  5.5  under  NUMBER  and  the  newly  obtained  LIGHT 
percent  T value  under  percent  T.  This  procedure  is  continued  until  values  have 
been  added  to  the  table  for  NUMBERS  4.5,  3.5,  and  2.5.  Additional  starting 
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points  such  as  7.8,  7.6,  etc.,  were  chosen  and  the  above  example  was  followed  to 
obtain  numerous  data  points  in  order  to  precisely  draw  the  emulsion  calibration 
curve.  Typical  data  points  used  to  establish  the  preliminary  and  emulsion  cali- 
bration curves  are  shown  in  Tables  1 and  2. 

The  emulsion  calibration  curve  just  completed  was  expanded  to  15"  x 10" 
semilog  paper.  The  ordinate  remained  the  same  scale;  however,  the  abscissa  was 
expanded  in  such  manner  that  two  curves  were  plotted,  one  over  the  other.  The 
NUMBER  or  abscissa  axis  scale  is  enlarged  so  that  the  NUMBERS  5 to  8 on  the  small 
calibration  curve  were  expanded  to  full  scale  for  the  new  curve.  Starting  with 
5,  each  major  division  is  in  tenths  and  each  minor  division  is  equal  to  0.02. 

The  top  abscissa  of  the  new  emulsion  calibration  curve  expands  the  scale  of  the 
small  calibration  curve  region  NUMBERS  2 to  5 to  full  scale.  The  enlarged  curve 
is  completed  using  values  from  the  table  used  to  construct  the  small  calibration 
curve . 

The  last  stage  of  the  calibration  process  involves  construction  of  a rule. 

A cardboard  approximately  28"  x 1.5"  is  used  as  the  starting  material  and  a line 
is  drawn  dividing  the  cardboard  in  half  in  the  long  direction.  Abscissa  axis 
NUMBERS  2.0,  3.0,  4.0,  5.0,  6.0,  and  7.0  are  superimposed  on  the  rule  from  the 
expanded  calibration  curve  and  marked  off  from  the  underside  of  the  long  line  of 
the  rule.  This  is  done  in  such  manner  that  approximately  abscissa  axis  NUMBER 
7.5  would  be  the  last  point  at  the  right  hand  end  of  the  rule. 


With  the  rule  aligned  so  that  NUMBERS  7,  6,  and  5 on  the  graph  paper  are  in 
line  with  their  corresponding  NUMBERS  on  the  rule,  a right  triangle  is  used  to 
help  superimpose  the  percent  T ordinate  values  of  the  large  emulsion  calibration 
curve  onto  the  rule  and  to  mark  off  percent  T readings  on  the  upperside  of  the 
straight  line  of  the  rule.  For  example,  using  2 percent  T,  one  finds  the  inter- 
section of  the  ordinate  point  2%  with  the  curve,  aligns  the  long  straight  edge 
of  the  triangle  at  this  point  while  the  base  of  the  triangle  is  parallel  with  the 
straight  line  of  the  rule,  an  alignment  mark  is  made,  and  the  NUMBER  2.0  is 
labeled.  In  this  manner,  percent  transmissions  between  2.0  and  approximately 
60.0  were  marked  off.  The  rule  was  then  shifted  so  that  abscissa  axis  NUMBERS 
2.0,  3.0,  and  4.0  on  the  calibration  curve  were  aligned  with  the  corresponding 
numbers  on  the  rule  and  the  remainder  of  the  percent  T scale  was  superimposed. 

The  rule  has  the  approximate  range  92.0  to  2.0  percent  transmission. 


Table  2.  TOPICAL  LOG  PERCENT  TRANSMISSION  AND 

Table  1.  TYPICAL  RAW  DARK  AND  LIGHT  NUMBER  VALUES  FOR  ESTABLISHING  THE  EMULSION 

DATA  POINTS  USED  TO  ESTABLISH  THE  CALIBRATION  EQUATION 


PRELIMINARY  CURVE  EQUATION 


Point 

Dark 

Light 

Point 

Log  Percent 
Transmission 

Percent 

Transmission 

1 

8.7 

24.5 

7 

0.4047 

2.54 

7 

5.1 

14.7 

6.6 

0.5671 

3.69 

14 

42.1 

77.1 

6.2 

0.7295 

5.36 

15 

65.1 

91.9 

5.8 

0.8917 

7.79 

19 

11.1 

29.1 

5.4 

1.0528 

11.29 

20 

59.5 

85.8 

5.2 

1 .1328 

13.58 

30 

49.7 

81 .1 

4.8 

1.2907 

19.53 

34 

30.5 

58.1 

4.3 

1 .4806 

30.  ?4 

53 

38.1 

73.3 

3.7 

1 .6859 

48.52 

73 

2.9 

6.2 

3.3 

1.7974 

62. 72 

2.5 

1 .9318 

85.4 
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After  establishing  the  rule,  one  is  then  ready  to  interpret  densitometer 
readings  of  lines  on  photographic  plates  in  the  calibrated  wavelength  region, 
corresponding  to  elements  under  study.  A percent  T for  the  known  and  for  the 
internal  standard  is  read  from  the  glass  emulsion.  Going  to  the  rule,  and  using 
an  ordinary  ruler  graduated  in  24  units  to  the  inch  which  corresponds  to  the 
linear  axis  spacing  on  the  semilogarithmic  paper  recommended,  the  relative 
intensity  difference  can  be  obtained.  By  setting  the  zero  on  the  left  end  of 
the  ruler  on  the  higher  percent  T,  one  then  notes  the  point  on  the  ruler  which 
corresponds  to  the  other  percent  T.  That  number  is  the  relative  intensity  for 
that  pair  of  percent  T values.  Continuing  this  process  for  pairs  of  known  and 
internal  standard  lines,  an  analytical  curve  can  be  established  on  semi  log  paper 
with  the  log  axis  labeled  CONCENTRATION  and  the  linear  axis  RELATIVE  INTENSITY. 
From  this  curve  unknown  concentrations  with  known  relative  intensities  can  be 
determined.  Typical  preliminary  and  emulsion  calibration  curves  and  a rule  are 
shown  in  Figures  1 to  3. 


HEWLETT-PACKARD  CALCULATOR  METHOD 

The  equation  of  the  preliminary  curve  can  be  determined  and  the  curve  gen- 
erated by  a plotter  by  the  use  of  the  polynomial  regression  procedure  outlined 
in  the  MATH  PAC  book  furnished  with  the  H.P.  computer.11  Following  this  method 
one  inputs  a set  of  data  points  (Figure  4),  and  coefficients  of  a polynomial  up 
to  the  20th  degree  using  a least-square  fit  may  be  produced.  The  DARK  percent  T 
and  LIGHT  percent  T data  used  to  draw  the  preliminary  curve  (see  method  done  by 
hand)  was  entered  into  the  H.P.  program.  Various  degrees  of  the  polynomial  were 


Figure  1.  Preliminary  curve, 
exemplifying  method  of 
successive  determination  of 
DARK  %T  and  LIGHT  %T 
values. 


I ).  98tO  A MATH  PAC.  Model  30,  Hewlett-Packard  Calculator,  v.  1. 


Number 


Condensed  rule 


examined  and  it  was  found  that  a 6th  degi’ee  gave  a sufficiently  high  correlation 
coefficient  for  use.  The  following  equation  was  used: 


where  Ti  = percent  transmission  of  light  line 
Tj  = percent  transmission  of  dark  line 
A c,  Ai ....  coefficients  of  polynomial  equation 
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H.  P.  CALCULATOR  OPERATIONS  CHART 


Tj  = Percent  transmission  of  light  line 
Td  = Percent  transmission  of  dark  line 


Figure  4.  H.P.  calculator  operations  chart. 
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The  data  outputted  contained  DARK  percent  T values  beginning  with  the  number 
1 and  incremented  in  units  of  1 to  80  along  with  corresponding  LIGHT  percent  T 
values  which  are  the  calculated  value  for  each  DARK  percent  T point.  Plotting 
the  data,  a smooth  curve  was  obtained.  This  curve  was  superimposed  on  a hand- 
drawn  curve  from  the  original  data  and  was  found  to  be  essentially  identical. 

Using  the  coefficients  of  the  6th  degree  polynomial,  one  generates  a table 
labeled  LIGHT  percent  T and  NUMBER  which  is  used  to  formulate  the  emulsion  cali- 
bration curve.  Referring  to  the  above  equation,  let  Td  equal  1,  solve  the  above 
equation  for  Tj  and  store  this  value  in  the  table  along  with  the  NUMBER  7.  The 
computer  sets  Td  = Tj  and  resolves  the  equation.  In  this  manner,  values  in  the 
table  for  Tj  are  established  for  NUMBERS  7,  6,  5,  4,  3,  and  2. 
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A linear  curve  is  assumed  between  the  data  points  NUMBER  = 7,  percent  T = 
value  found,  and  NUMBER  = 6,  percent  T = value  found.  The  straight  line  region 
is  extended  to  percent  T = 1,  NUMBER  = calculated  value  from  the  equation  of  the 
straight  line.  Solving  this  straight  line  equation,  the  computer  is  then  able 
to  take  NUMBER  values  along  this  line  and  obtain  corresponding  Tj  values.  For 
example,  using  the  NUMBER  value  7.5,  a corresponding  value  is  obtained  from 
the  equation  of  the  straight  line.  This  is  set  equal  to  Td  and  placed  in  the 
previous  polynomial  equation: 

Ti  = Aq  + Ai  (Td)  +....A6  (Td)6 

and  a corresponding  Ti  value  is  received.  In  the  table  T^  and  NUMBER,  the  T^ 
value  obtained  and  the  NUMBER  6.5  (one  less  than  7.5),  are  stored.  In  similar 
fashion,  NUMBER  values  from  6.9  to  2.1  in  increments  of  one  tenth  are  obtained. 

The  data  from  the  table  LIGHT  percent  T versus  NUMBER  is  then  used  to  gen- 
erate the  coefficients  for  the  equation: 

NUMBER  = B0  + Bi  (log  TJ  + B6((log  Tx)6. 

A 6th  degree  polynomial  furnished  a high  correlation  coefficient. 

Using  the  above  equation,  percent  T for  a known  element  (obtained  from  the 
density  of  the  line  on  the  glass  plate  emulsion)  is  substituted  for  the  Ti  value. 
Solving  the  equation  a value  for  NUMBER  is  obtained.  The  procedure  is  repeated 
for  a percent  T for  an  internal  standard  and  another  NUMBER  is  obtained.  The 
computer  subtracts  these  two  NUMBERS  and  their  difference  is  called  RELATIVE 
INTENSITY.  As  outlined  in  the  method  done  by  hand,  this  process  is  continued 
for  pairs  of  known  and  internal  standard  lines.  An  analytical  curve  is  estab- 
lished on  semilog  paper  with  the  log  axis  labeled  CONCENTRATION  and  the  linear 
axis  RELATIVE  INTENSITY.  From  this  curve  unknown  concentrations  with  known 
relative  intensities  can  be  determined.  Program  used  is  presented  in  Appendix  B. 


UN  I VAC  1108  COMPUTER  METHOD 

Appendix  C lists  the  UNIVAC  1108  program.  The  only  input  necessary  for  the 
program  is  the  initial  DARK/LIGHT  data  read  from  the  emulsion-covered  glass  plate 
used  for  the  calibration.  Using  these  values,  the  computer  can  calculate  coef- 
ficients of  a polynomial  equation  to  any  practical  degree  required,  and  from  the 
values  obtained,  the  lowest  degree  polynomial  with  good  correlation  coefficient 
was  chosen.  The  program  then  outputted  least-square  data  of  DARK/ LIGHT  points, 
that  is  DARK,  LIGHT,  and  LIGHT  FITTED,  from  which  the  preliminary  curve  may  be 
drawn.  If  desired,  a plot  of  this  curve  can  be  obtained. 

Following  the  procedure  outlined  in  the  Hewlett-Packard  Calculator  Method, 
a table  labeled  LIGHT  percent  T and  NUMBER  is  generated  to  formulate  the  coeffi- 
cients of  the  emulsion  calibration  equation  used  as  the  equation  of  the  emulsion 
calibration  curve. 

NUMBER  = B0  + Bx  (log  TJ  + Bn  (log  Ti)n 


7 


1 

r 


where  T*  = percent  transmission  of  light  line 

B0....Bn  = coefficients  of  polynomial  equation. 

After  determining  the  coefficients  of  the  above  equation  with  the  UNIVAC  1108, 
the  equation  can  be  used  with  many  small  desk  top  calculators  that  do  not  have 
tie  memory  needed  to  solve  a 6th  degree  polynomial.  By  inserting  the  log  of  the 
percent  transmission  of  an  elemental  line,  one  gets  the  value  for  NUMBER. 

Assuming  a calculator  is  not  available,  using  the  above  equation,  a chart 
of  percent  transmission  and  NUMBER  may  be  generated.  Percent  T was  varied  in 
increments  of  0.05  from  2 to  20  and  in  units  of  0.1  from  20  to  90.  The  left 
vertical  column  of  the  chart  lists  changes  in  percent  transmission  by  one  unit, 
far  example,  from  20  to  21.  The  top  row  of  the  chart  varies  percent  T by  either 
0.05  or  0.1  depending  on  the  percent  transmission  range.  The  body  of  the  chart 
contains  NUMBER  values.  For  example,  if  an  intensity  of  an  element  of  a standard 
sample  line  was  measured  at  50.5%  T and  that  of  an  element  internal  standard  line 
at  20.6%  T,  one  would  find  the  corresponding  NUMBER  value  by  proceeding  horizon- 
tally along  the  chart  from  the  extreme  left  column  containing  the  number  50  as 
the  first  integer  until  one  reaches  the  column  headed  by  0.5.  This  NUMBER  value 
is  then  recorded.  The  procedure  is  repeated  for  the  other  transmission  value  of 
20.6.  The  difference  in  NUMBER  values  is  the  relative  intensity. 

Vari  .s  harts  were  attempted,  including  one  in  which  two  percent  trans- 
missions v.foi  jxample,  from  the  unknown  element  line  and  that  of  the  internal 
standard)  could  be  found  in  the  first  column  and  first  row  of  the  chart.  The 
point  of  intersection,  going  horizontally  and  vertically  along  the  chart  from 
the  two  percent  T values,  is  the  relative  intensity  and  the  need  of  taking  the 
difference  in  NUMBER  for  each  percent  T pair  is  eliminated.  To  accomplish  this, 
however,  the  chart  became  very  voluminous  and  it  was  found  that  the  previously 
described  chart  was  easier  to  use. 


CONCLUSION 

A computer  program  that  follows  closely  the  ASTM  graphic  method  has  been 
developed  for  use  on  a Hewlett-Packard  9830-A  and  a UNIVAC  1108  time-sharing 
computer.  For  laboratories  that  do  not  have  continuous  access  to  these  computers, 
one  can  obtain  the  coefficients  to  the  equation  relating  percent  transmission  to 
relative  intensity  and  then  use  the  equation  on  small  memory  desk  top  calculators 
to  find  relative  intensities. 

The  program  used  in  the  UNIVAC  1108  will  output  a table  of  percent  trans- 
nission  versus  NUMBER  which  is  used  to  obtain  the  relative  intensity  for  a pair 
cf  percent  transmission  values. 
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APPENDIX  A.  SPECTROGRAPH  1C  PARAMETERS 


Equipment  0 

Baird  3-meter  spectrograph,  2.75  A/mm  reciprocal  linear  dispersion 
Jarrel  1 -Ash  1/arisource 
Jarrell -Ash  Microphotometer 

Instrumental  Parameters 

Spark-18,000  V.  A.  C.,  3 breaks  per  half-cycle,  6 to  6.5  amperes 

Analytical  gap-3  mm 

Auxi 1 iary  gap-6  mm 

Sample  electrode-Ul tra  Carbon  106 

Counter  electrode-Ultra  Carbon  105  U 

Pre-spark  duration-20  sec. 

Photographic  Parameters 

Plates-Eastman  Kodak  SA  1 
Developer-Eastman  Kodak  D 1 9 
Developing  time-  3.5  min. 

Developimg  temperature  20°C 

0 

Spectral  range  of  calibration  2750-3A50  A 
Spectral  lines  read  in  range  of  3000-3200  A 
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APPENDIX  B.  HEWLETT-PACKARD  CALCULATOR  PROGRAM  FOR  GLASS  PLATE  EMULSION  CALIBRATION 


10 

For  1-1  to  22 

30 

End 

20 

C [1]  - B ll]  = 0 

A0 

Disp  " Wrong  No.,  Dark,  Light 

30 

Next  1 

50 

Input  1,  Y,  B [2] 

AO 

For  1 = 23  to  253 

60 

G [1]  = Y 

50 

C [1]  = 0 

70 

H [1]  - B [2] 

60 

Next  1 

80 

Print  "Delete:  Dark  = "G[|]M  Light 

70 

B [1]  - 1 

90 

Y = FNX  (-1) 

80 

W - N = S 1 = S2  = S3  = SA  = S5  = II  =0 

100 

D i sp 

90 

Disp  " Max  Degree  Polynomial 

110 

End 

100 

Input  D2 

f3 

1 10 

Print  " Max  Degree  Polynomial  = " 02 

10 

11=0 

120 

D 1 = D2 

20 

Print  " Polynomial  Degree  = "Dl 

130 

Print 

30 

Print  "Relative  1 Log  T Light  " 

1A0 

End 

A0 

II 

CNJ 

X 

50 

X = 1 

10 

Format  F A.0,  2 F 12. A 

60 

XI  = 8 

20 

Print  "Enter  Prelim  Calibrn  Data" 

70 

If  Xl<  = 2 then  210 

30 

If  N * 0 then  70 

80 

Y = B [D 1 +1] 

AO 

Print  "No."  Tab  10"  Dark  "Tab  22  "Light" 

90 

For  J * Dl  to  1 Step  -1 

50 

1 1 - 0 

100 

Y = Y * X + B[j] 

60 

11=11+1 

no 

Next  J 

70 

Disp  "Dark,  Light 

120 

X-Y 

80 

Input  Gill],  H [ 1 1 ] 

130 

If  X> 1 00  then  210 

90 

Y * H [II] 

1 AO 

XI  = Xl-l 

100 

B [2]  = G[ll] 

150 

If  Xl<  = 2 then  210 

110 

Write  (15,  10)  N+l , G [II],  H [ll] 

160 

11=11+1 

120 

Y - FN  X 1 

170 

H [II]  = XI 

I3C 

Go  to  60 

180 

G [11]  • L G T I 

IA0 

End 

190 

Print  XI,  G [II],  X 

200 

Go' to  70 

10 

If  W = 0 then 

210 

X2  = X2  - 0. 1 

20 

Disp  " Not  a 1 lowed" 

220 

If  X2  = 0 then  290 

230 

11=11+1 

"H  [I] 


'JSl- 


10 


-M -fy 


370 

D2  = D2  + 1 

710 

Next  J 

t 1 

380 

For  J = 1 to  D2 

720 

Next  1 

1 f 

s 

390 

Print  C [P] 

730 

c [1]  - 1/  c [1] 

it  00 

C [P]  - SQR  C [P] 

740 

Go  to  60 

it  1 0 

For  1 = 1 to  D2  - J+l 

f7 

it20 

C (P+l]  = C [P+l]  7 C [P] 

10 

e [22]  = Dt 

430 

Next  1 

20 

Print  " Store  X = F ( Loo 

440 

R = P + 1 

30 

Di sp  " Store  Fi le  0 = " ; 

l ! 

450 

S - R 

40 

Input  1 

t 

460 

For  t = 1 TO  02  - J 

50 

Pr int  " Store  File  V " 1 

470 

P = P + 1 

60 

Store  Data  1 , B 

480 

For  M = 1 to  D2  ♦ 2 - J - L 

70 

End 

490 

C [ R+M  - 1]  + C [ R+M-l]  - C [P]  * C [P+M-l] 

f8 

500 

Next  M 

10 

Print  " Oiff  X Evaluation 

510 

R - R + M - 1 

20 

Disp  " T (STD) , T « "J 

520 

Next  l 

30 

1 nput  Tl , T2 

1 

530 

P = S 

40 

Print  " T (STD)  = "Tl 

540 

Next  J 

50 

Print  " T - " T 2 

550 

T - ( D2  + 1 ) ' (02  + 2)  /2 

60 

Tl  = LGTTI 

560 

For  1 = 1 to  02  - 1 

70 

T2  = LGTT2 

570 

T = T - 1 - 1 

80 

XI  - X2  - B [Dl  + 1] 

580 

C [T]  " 1 / C [T] 

90 

For  J = Dl  to  1 Step  - 1 

[ 

590 

For  J = 1 to  02  - 1 

10C 

XI  * XI  * Tl  + B [J] 

600 

P * 02  + 1 - 1 - J 

110 

X2  - X2  * T2  + B (J] 

■ 

610 

P « P * (02+1  - (P-l)  /2)  - 1 

120 

Next  J 

- 

i> 

620 

R - P-J 

130 

Print  " X (STD)  = " xl 

! 630 

S = 0 

140 

Print  " X « " X2 

«■ 

V, 

640 

II  = 1 + J + 1 

150 

Print  " X (STD)  - X « "XI 

• 

650 

Y * P 

160 

Go  to  20 

660 

For  K * 1 to  J 

170 

End 

670 

Y * Y f U - K 

c 

'3 

•• 

680 

S = S - C [R+K]  * c [Y] 

J 10 

Disp  "Re-plot  Axis"; 

■ 

690 

Next  K 

20 

Input  Z9 

700 

C [P)  - S / C [R] 

30 

If  29  = 0 then  390 

, 

V.  * i' 

t. 

«|v 

. « • 

12 
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40 

Disp  "X  (Max)  = 

390 

For  X = A2  to  B2  Step  C2/10 

50 

Input  B2 

400 

Y • B [01  ♦ 1] 

60 

Oisp  "X  (MJn)  = 

410 

For  J = 01  to  1 Step  -1 

70 

Input  A2 

420 

Y = Y * X + B [J] 

80 

Oisp  "Y (Max)  * 

430 

Next  J 

90 

Input  Bt 

A40 

Plot  X,Y 

ion 

Disp  "Y  (Min)  * " ; 

450 

Next  X 

110 

Input  A1 

460 

Disp  "Plot  Data  Points"; 

120 

C2  - ( B 2 - A2)  / 10 

470 

Input  29 

130 

Cl  - (81  - Al)  / 10 

480 

If  29  = 0 then  550 

HO 

Scale  A2,  B2  + 2*C2,  Al,  El  , + 2*C1 

490 

PEN 

150 

Plot  A2  + (B2  - A2)  12.2,  Al  + Cl  /4 

500 

For  K = 1 to  1 1 

160 

Label  (*,  3,  2,  0,  2/3) 

510 

Plot  G [K] , H [K] 

170 

Letter 

520 

CPLOT  -0.3,  -0.3 

180 

Plot  A2  + C2/4,  Al  + ( B 1 -Al )/2.2 

530 

Label  (*)  "0" 

190 

Label  (*,  3,  2,  PI/  2,  2/3 

540 

Next  K 

200 

Letter 

550 

PEN 

220 

Offset  C2,  Cl 

560 

Plot  -1.2  * C2 , -Cl 

230 

XAXIS  Al,  C2,  A2 , B2 

570 

PEN 

2ltO 

YAXIS  B2 , Cl , Al , B1 

580 

Go  to  10 

250 

YAXIS  A2,  Cl , Al , B 1 

f,l 

260 

XAXIS  B1 , C2,  A2,  B2 

10 

Format  2F  12.4 

270 

Label  (*,  2,  1.7,  0,  2/3) 

20 

Print  Tab  9 "X  Dark  "Tab  21' 

280 

For  Y = Al  to  B1  Step  2*0 

30 

For  1 * 1 to  79  Step  2 

290 

Plot  A2 , Y,  1 

40 

Y « B [01  + 1] 

300 

CPLOT  -5.A,  -0.34 

50 

For  J • 01  to  1 Step  -1 

310 

Label  (380,  2,  1.7,  0,  2/3)  Y 

60 

Y - Y*  1 + B [J] 

320 

Next  Y 

70 

Next  J 

330 

For  X = A2  to  B2  step  2*C2 

80 

Write  (15,  10)  1,  Y 

3W 

Plot  X,  Al,  1 

90 

Next  1 

350 

CPLOT  -3.8,  -1 

100 

For  1 = 80  to  100  Step  1 

360 

Label  (380,  2,  1.7,  0,  2/3)  X 

no 

Y - B [Di  + 1) 

370 

Next  X 

120 

For  J = D 1 to  1 Step  -1 

380 

Format  F 5*1 

130 

Y - Y * 1 + B [ J ] 

**'•••>  VI 


APPENDIX  C.  UNIVAC  1108  PROGRAM 


1 

105 

103 
1 
2 

104 
1 


1 

10 

2 

101 


1 


82 

81 

15 

75 

80 


DIMENSION  X(1200),  Y(1200),  N(200) , Tl(200),  T2(200),  T3(200), 

C(ll),  ALPHA (10)  , BETA (10)  , 

DIMENSION  S(4),  YF(400),  A(14) 

FORMAT  ( III  29H  OUTPUT  FROM  SUBROUTINE  FITY.//7H  IND  3 = I 3///3X 
2H  IX,  18X  2HY,  14X  9H  Y FITTED  //  (1H  1P3E  20.7)) 

FORMAT  ( / / / 3 1 H OUTPUT  FROM  SUBROUTINE  ORTHLS.//7H  IND  1 = 13//2X 
2HI,  13X  2H  C,  16X  6H  ALPHA,  14X  5H  BETA/1H  14,  1PE  20.7/ 

(IK  14,  1P3E  20.7) j 

FORMAT  (///  3 OH  OUTPUT  FROM  SUBROUTINE  COEFS.//7H  IND  2 = I3//13X 
2HI , 17X  SH  A(I)//(1H  114,  1PE  26.7)) 

P = 3.141592654 
READ  1,N,K 

READ  1,  NUM,  NUMA,  DUM 
FORMAT  (IS,  15,  F10.0) 

DO  10  I = 1,  N 
READ  2,  X(I),  Y(I) 

CONTINUE 

FORMAT  (F10.0,  F10.0) 

L = 0 
J = 0 

CALL  ORTHLS  (X,  Y,  W,  N,  L,  J,  C,  ALPHA,  BETA,  K,  T1 , T2,  T3,  IND1) 

IT  * 0 

WRITE  (6,  103)  IND  1,  II,  C(l),  (II,  C (I I + 1),  ALPHA  (II),  BETA  (II), 
II  * 1,  K) 

DO  15  MM  = 1,K 
KF  = MM 

CALL  FITY  (X,  N,  J,  C.  ALPHA,  BETA,  KF,  YF,  T] , T2,  IND  3) 

WRITE  (6,  105)  IND  3,  (X ( II  . Y(I),  YF ( I ) , I = 1,  N) 

CALL  COEFS  (J,  C,  ALPHA,  BETA,  KF,  A,  T1 , T3,  IND  2) 

IT  * 0 

WRITE  (6,  104)  IND2 , II,  A(l),  II,  A(II  ♦ 1), 

II  - 1,  KF) 

FORMAT  (1H1) 

FORMAT  (1H  3E  15.6) 

MB  = NW  +1 
CONTINUE 
PRINT  82 

DO  80  IA  = 1,  N 
ANS  = A(l) 

DO  75  II  = 2,  MB 

FIB  = II  -1 

ANS  * ANS  + A(II)  * X f I A) * FIB 
IF  (IT.  EQ.3)  X ( I A)  = EXP  (X( IA) ) 

PRINT  81,  X(IA) , Y (IA) , ANS 
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nm 
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IF,  (IT.F.Q.3)  GO  TO  102 
PRINT  82 
AB  = 1.0 

SAVE  = SOLVE  (A,  AB,  MB) 

PRINT  85,  SAVE,  AB 

SAVE  = SOLVE  (A,  SAVE,  MB) 

PRINT  85,  SAVE  1,  SAVE 
CON  = ALOG  (SAVE) 

CON  1 = ALOG  (SAVE  1) 

D » CONI -CON 

Yf 1)  = ALOG  (SAVE) 

X(l)  * 1 

S(l)  = SAVE 

DD  * D/10 

DO  6S0  I * 2,  10 

SCI)  = EXP  (CON  + DD  * (I  - 1)) 

X(I)  = 7.  -(I-l)*.l 

650  Y(I)  = ALOG  (SCI) J 
L = 10 

DO  651  J * 1,  10 
AB  = S(J) 

DO  652  JJ  = 1,4 
L = L + 1 

ANS  = SOLVE  (A,  AB,  MB) 

Y(L)  = ALOG  (ANS) 

X(L)  = X(J)  -1.0  * JJ 
IF  (ANS.  GE.  100.)  GO  TO  653 
PRINT  85,  ANS,  AB,  X(L) 

652  AB  = ANS 
GO  TO  651 

653  L ■ L-l 

651  CONTINUE 
IT  = 3 

N = L 

DO  107  J = 1,  N 
TEM  = X(J) 

X(J)  = Y (J) 

107  Y (J)  = TEM 

GO  TO  101 
102  CONTINUE 

NUMB  = NUM  + 10 


NUMC  = NUMA  + 1 
DO  111  1 = 1,  NUMB 
X(I)  = DUM 

IF  (I.  GT.  NUMC)  X(t)  = X(l)  + .05 


V 


I 
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r ■ 


■ I 


AB  - ALOG  (X( I ) ) 

Y(I)  = SOLVE  (A,  AB,  MB) 

111  DUM  = X(I)  * .05 
PRINT  17 
WRITE  (6,  18) 

DO  10  I ■ 1,  NUMA,  20 
JJ  ==  I ♦ 19 

40  WRITE  (6.  19)  X C I ) , (Y(J),  J > I,  JJ) 
NUMB  = JJ  * 1 
I PRINT  = NUMB  ♦ 500 
PRINT  17 
WRITE  (6,  13) 

DO  ISO  I = NUMB,  NUM,  10 
JJ  = I ♦ 9 

IF  Cl-  GT.  IPRI.NT)  GO  TO  16 
GO  TO  150 

16  IPRI.NT  = I * 500 
PRINT  17 
PRINT  13 


150 

WRITE 

(6,11)  X(I), 

(Y (J)  , 

J = 

I, 

JJ) 

19 

FORMAT 

(1H  21F6 

.2) 

18 

FORMAT 

(16  X, 

,05' 

3X,  • 

.10' 

3X , 

'.15' 

3X, 

' .20  3X,  ' , 

.25' 

3X , '. 

30', 

1 

3X,  ' . 

35'  3X , ', 

.40' 

3X , • 

.45' 

3X , 

'.50' 

' 3X, 

'.55'  3.X , 

'.60'  3X,  '.65 

2 

3X, 

70'  3X , ' . 

.75' 

3X , 1 

.80' 

3X , 

’.85' 

' 3X, 

'.90'  3X , 

1 .95' 

11 

FORMAT 

(1H  F6.2, 

, 10F8.2) 

13 

FORMAT 

(1H  19X 

M' 

1 , 6X  ' 

.2', 

6X 

’ .3' , 

6X  ' 

.4’,  6X  '.5 

',  6X  '.6', 

6X  '.7',  6X  ’.8',  6X  '.9') 

17  FORMAT  (1H1) 

85  FORMAT  (1H  'Y  1 = 1 F13.8,  5X,  'XI  = ' F13.8,  E14.8) 

STOP 

F.ND 

LEGEND 

N = No.  OF  EXPERIMENTAL  POINTS 

K = MAXIMUM  ORDER  OF  EQUATION  FOR  POLYNOMIAL  LEAST  SQUARES  FIT. 

NUM  = TOTAL  NUMBER  OF  PERCENT  TRANSMISSION  vs  RELATIVE  INTENSITY  VALUES 
IN  TABLE 

NUMA  = TOTAL  NUMBER  OF  PERCENT  TRANSMISSION  VALUES  OF  INCREMENT  .05  IN  TABLE 

DUM  = INITIAL  VALUE  OF  PERCENT  TRANSMISSION  PRINTED  IN  TABLE  (2.0  USED  IN 
THIS  WORK) 

X(I)  = EXPERIMENTAL  DARK  PERCENT  TRANSMISSION 
Y(I)  = EXPERIMENTAL  LIGHT  PERCENT  TRANSMISSION 


17 


i 


No.  of 
Copies 


To 


No.  of 
Copies 


To 


* 
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2 Metals  and  Ceramics  Information  Center,  Battel le  Columbus 

Laboratories , 505  King  Avenue,  Columbus,  Ohio  43201  1 

1 

12  Commander,  Defense  Documentation  Center,  Cameron  Station,  1 

Building  5,  5010  Duke  Street,  Alexandria,  Virqinia  22314  1 

1 

Office  of  Chief  of  Research  and  Development,  Department  of 
the  Army,  Washington,  D.  C.  20310 
1 ATTN:  CRDPES 

1 

I ( ommander , Army  Research  Office,  P.  0.  Box  12211,  Research 
Triangle  Park,  North  Carolina  27709 

1 


Commander,  U.  S Army  Material  Development  and  Readiness  1 

Command,  5001  Eisenhower  Avenue,  Alexandria,  Virqinia  22313  ? 

1 ATTN:  DRCQA  2 

2 DRCQA-E  1 

2 DRCQA- P I 

1 DRCLDC,  Mr.  R.  Zentner  1 

1 DRCRL  1 

1 DREDE-EA  2 

1 DRCDE-E  1 

1 DRCRP-0IP 

1 

Commander,  U.  S.  Army  Electronics  Command,  Fort  Monmouth, 

New  Jersey  07703  2 

2 ATTN:  DRSFL-CB 

1 DRSEL-PP-PO 

1 DRSEL-WM  1 

1 DRSEL-RD-GT  1 

1 DRSEL-PA-C  1 

1 

' ommander,  U.  S.  Army  Missile  Command,  Redstone  Arsenal,  1 

Alabama  368^9  1 
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1 DRSMI-RSM,  Mr.  E.  J.  Wheel ahan 

1 DRSMT-RTR,  Mr.  H.  T.  Lawson  ? 

2 DRSMI-Q  1 

1 DRSMI-M  2 

Commander,  U.  S.  Army  Mobility  Equipment  Research  and 
Development  Command,  4300  Goodfellow  Boulevard,  St.  Louis, 

Missouri  63120  1 

1 ATTN:  DRSMF-Q 

1 0RSME-QP 

1 DRSME-QR 

1 DPSME-QE  1 

1 DRSME-M 

1 DRSME-P 

? DRSME-R 

2 

Commander,  U.  S.  Army  Armament  Command,  Rock  Island, 

Illinois  61201  2 

1 ATTN:  DRSAR-SC , Dr.  C.  M.  Hudson 

1 DRSAP-PPW-PB,  Mr.  Francis  X.  Walter  1 

2 Technical  Library  1 

1 

2 Commander,  U.  S.  Army  Natick  Research  and  DevelODment  1 

Command,  Kansas  Street,  Natick,  Massachusetts  01760  1 

i ATTN:  PRXNM-QE  1 

1 DRXNM-GE  1 

1 
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Development  Command,  28251  Van  Dyke  Avenue,  Warren,  1 

Michigan  48090  1 

1 ATTN:  DRDTA-PXMT,  Mr.  J.  Dudzinski  1 
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2 
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Proving  Ground,  Maryland  21005  1 

2 ATTN-  DPSTE-TA-A 

Commander,  Aberdeen  Proving  Ground,  Maryland  21005  1 

1 ATTN;  STFAP-MT,  Mr.  J.  M.  McKinley 

1 steap-tl 

1 


Commander,  Edqewood  Arsenal,  Edgewood,  Maryland  21010 
ATTN:  SAREA-TSP 

SAREA-QAE 
SAREA-QAI 
SAREA-QAP 
SAREA-QAIP 

Commander,  U.  S.  Army  Foreign  Science  and  Technology  Center, 

220  7th  Street,  N.  E.,  Charlottesville,  Virginia  22901 
ATTN:  DRXST-SD2 

Commander,  Frankford  Arsenal,  Philadelphia,  Pennsylvania  19137 
ATTN:  SARFA-P3300 
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SARFA-Q2100 
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Adelphi,  Maryland  20783 
ATTN:  DRXDO-EDE 

Commander,  Picatinny  Arsenal,  Dover,  New  Jersey  07801 
ATTN:  SARPA-RT-S 

SARPA-VA6 , Mr.  H.  DeFazio 
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SARRI-LEQ 
SARRI-LRM 
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